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Abstract. The galaxy cluster A3667 was observed using the Two-degree Field (2dF) multifibre spectroscopic 
system on the Anglo-Australian Telescope in a program designed to examine the velocity structure in the region. 
Specifically, we sought evidence from the optical data for the putative cluster merger believed to be responsible 
for the observed radio and X-ray emission. We present 184 new redshifts in the region, of which 143 correspond 
to member galaxies of A3667. We find the cluster velocity distribution to be well modelled by a single Gaussian 
in agreement with previous results. In addition, new redshift-selected isodensity plots significantly reduce the 
prominence of the previously reported subgroup to the north-west of the main cluster. Instead, we find the 
galaxy distribution to be elongated and well mixed, with a high velocity dispersion and no significant evidence 
for substructure. These results are consistent with the axis of the proposed merger being close to the plane of the 
sky. 
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^ , 1. Introduction 



lected by appare nt magnitude, an d a similar conclusion 



The AGO galaxy cluster, A3667 (|Abell et al.lll989l ). is a 
rich. X-ray-luminous cluster, which is arguably one of the 
most interesting and, consequently, most studied clusters 
in the southern sky. It has the distinction of being one 
of only a handful o f galaxy clusters to h ave a cold gas 
front at the centre (jVikhlinin et al.l 120011 ) and the only 
cluster yet seen with two bright , diffuse, radio emission 



was reached by Sodre et al. I (|l992h . 



ROSAT observations revealed a high X-r ay luminos- 
ity L x, consistent with the Lx — <J relation (jXue fc Wu 



20000. of 8.74 X 10^ "^ erg s^^ in the 0.4-2.4 keV range 



regions straddhng the X-ray gas ( Rottgering et al.|[r997,) 
Our physical picture of the cluster is becoming clearer as 
more mult i wavelength data are obtained, but many basic 
questions about the dynamical history of A3667 remain 
unanswered. 

A3667 is relatively close, with a published mean red- 
shift of 0.055 and a large velocity dispersion (ctv) of 

1000-1500 km s^^ depending on the number of galaxies 

I , ^ I ■ 

used and the redshift limits applied dG irardi et al. 199 _ 

Melnick fc Quintaiialll981 ; Proust et al.iil988.; .Sodre et al 



( Ebeling et al.|[l996h . making it one of the brightest X-ray 
sources in the southern skjo. Further, the X-ray isophotes 
are distorted in the direction of the reported bimodal opti- 
cal distribution and observations from the Einstein satel- 
lite are reported as sh owing clear evidence of substruc- 



ture llSodre et al.lll992l ). Observ ations of the temp e rature 



19920 . The Abell richness class is 2 and the Bautz- Morgan 



type is I-II (jAbell et al.lll989[). ind icating a high galaxy 



of the X-ray gas with ASCA by iMarkevitch et al.l (|l999t ) 
showed the central part of the cluster to be cool and the 
X-ray surface brightness profile along the elongation axis 
suggests a shocked region is present. Recent Chandra ob- 
servations confirmed this interpretation, revealing the cen- 
tral part of the cluster to contain a cold g as front moving 
through the warm intracluster medium ( Vikhlinin et al.l 
200 ih . 

Perhaps the most dramatic features of A3667 are 
the two extended, sy mmetrically located reg ions of dif- 



density in this region. [ Proust et al. | (11988,) presented ev- ^^^^ ^^^^^ emission (jRottgering et al.l [l997h . first de 



idence for a bimodal spatial distribution of galaxies se- ^^^^^^ -^^ §43 MHz images obtained with the Molonglo 
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Observatory Synthesis Telescope, and then confirmed and 
studied in detail at 1.4 and 2.4 GHz with the Australi a 
Telescope Compact Array, ATCA (| Johnston-Hollittll2003f ) . 
The combination of these features has been inter- 
preted as evidence that the cluster is observed in a 



post- m erger state after c l ose co r e passage (IKnopp et al. 



R5ttgering et al.lll997l : IVikhlinin et al.ll2001 ) 



1996t iMarkevitch et ail Il999l: iRoettiger et al.l Il999l: 



Given the strong indicator of merger activity present 
in the radio, X-ray and previous optical data we undertook 
a set of spectral observation of the objects in the cluster 
to further test for signatures of a cluster-wide merger in 
the optical population of A3667. 

In Section 2 we discuss the signatures of mergers that 
have been observed in galaxy clusters in the optical and 
present a detailed summary of the previous optical data 
collected for A3667. Sections 3 and 4 discuss the Two- 
degree Field (2dF) observations and data reduction. The 
major results are discussed in Section 5 and the discussion 
and conclusions are presented in Section 6. 

2. Previous Optical Observations of A3667 

Melnick fc Quintanal ( 1981 ) obtained redshi fts for 16 



member galaxies of A3667. Following this, iProust et al 



( 1988[ ). hereafter referred to as P88, obtained a further 31 
redshifts for galaxies in the field of A3667, which they com- 
bined with the previous results. Unfortunately, this com- 
bined sample has sin ce been found to c ontain some unre- 
liable redshift values (ISodre et al.lll992n . P88 also posited 



a bimodal spatial distribution of galaxies in the field of 
A3667 based on an isodensity plot of 423 galaxies with 
TTiB < 19. Two groupings can be seen in this plot. The 
main group was assigned a redshift of 0.0550, being that 
of the galaxy closest to the centre of A3667, while the sub- 
group to the north-west was suggested to be at z = 0.0563, 
based on the redshift of a nearby galaxy. It should be noted 
that nearly 90% of the galaxies in this analysis had un- 
known redshifts and there is an inevitable contribution 
fro m foreground and background objects. 



Sodre et al.l (|l992f ) expanded on the P88 data by mea- 



suring redshifts for 128 galaxies in the vicinity of A3667. 
This sample, hereafter referred to as S92, included reob- 
servation of all but five of the galaxies in the P88 sample 
and revealed some large discrepancies in redshift. As nei- 
ther study published spectra this has made it difficult to 
determine which was the more reliable; in general the S92 
sample has smaller quoted errors and is assumed to be 
more reliable. Other data in the literature found using 
a search with the NASA Extragalactic database (NED) 
have provided a further 28 redshifts for objects in the S92 
sample, bringing the total number of redshifts available in 
the region of A3667 to 161. Of these, 155 fall in the range 
0.0 44 < z < 0.06 8 and a re likely to be cluster members. 
Girardi et al. I (|l996h studied A3667 as part of a sam- 



ple of 37 galaxy clusters used to examine the relationship 
between velocity dispersion, CTv, and X-ray temperature, 
T. They argued that since both CTv and T were related to 



Table 1. Mean redshifts and velocity dispersions for 
A3667. 



Reference 



N Mean 



(kms 



Proust et al. (19881 
Sodre et al. (1992) 
Girardi et al. (1996) 
Girardi et al. (1998') 



This paper 



0.054±0.003 
0.0553±0.0004 

0.0566 



1606±288 
1197±79 
1208 



-1-95 
-84 

971+''^ 



-47 



231 0.0555±0.0003t 1102±46 



* From the table in lProust et al.1 (Il988l ). 

* The mean of all 231 available redshifts; the best-fit 
Gaussian to the data gives 1 = 0.0553 (see Section 6). 



a cluster's gravitational potential there would only be a 
strong correlation between these quantities for clusters in 
which both gas and galaxies were in dynamical equilibrium 
with the clu ster potential. Using the 123 redshifts from the 
S92 sample, iGirardi et al.l (|l996l ) examined A3667 for sub- 
structure using both a velocity distribution and velocity 
gradient analysis. They found neither significant multiple 
peaks in the velocity d istribution, nor a veloc i ty gra dient. 
Using the method of iDressler fc Shectman ( 19881 ) they 
computed the probability of substructure as 0.697 where 
a result above 0.990 was considered significant. Despite 
the lack of statistical evidence, their isodensity plot was 
reminiscent of the P88 plot and suggested some substruc- 
ture to the northwest of the cluster centre. The velocity 
dispersion was determined to be 120 8j^s4 km s~^. 

Expanding on their earlier work, Girardi et al.l ( 19981 ) 
used 154 of the then known 155 redshifts for A3667 to 
compute a maximum cluster radius of 2.22h~^ Mpc, a 
mean galactocentric redshift of 0.0 566 and a glo b al ve - 
locity dispersion of 971^47 km s^^. iGirardi et al. (1998) 
also obtained a virialised radius of 1.94/i~^ Mpc in which 
152 of the galaxies with previously measured redshifts 
are contained. The total virialized mass was calculated 
to be 15.98tl^jl x 10^'' M©. They further concluded that 
the velocity distrib ution was isotropic. In comparison 
Knopp et all (|l996h give the X-ray radius of the clus- 
ter to be Rx — 0.80h~^ Mpc and the X -ray mass at 



M 



X 



3.50: 



-0.30 
-0.40 



X 10^-* Mq 



Girardi et all |T998) obtain 
Q within this 



an optical mass of Mq = 4.89lo;f2 x lO" M, 
X-ray radius. Assuming that the cluster galaxies follow 
the gravitational potential, the good agreement between 
these mass estimates stro ngly suggests that th e cluster is 
in dynamical equilibrium ( Girardi et al. 1998f) . 

The mean redshifts and velocity dispersions obtained 
in the studies cited above are summarized in Table [T] 



3. 2dF Observations 

In order to increase the redshift sample and examine 
the distribution of galaxy types in the cluster, observa- 
tions were undertak en with the 2dF spectrographic system 
(iBailev et al.l l2002h on the AAT. SuperCOSMOS plate 
scans were used to obtain a list of suitable galaxies with 
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15.5 < me < 18.3. As there is a marked discrepancy be- 
tween the redshift values for some objects observed by P88 
and S92 we decided to include a subset of the S92 cata- 
logue in the target sample to check for any systematic dif- 
ferences in redshift. The opti cal counterparts of the k nown 
radio galaxies in the region (| Johns ton-HoUittl l2003h were 
also selected. This produced a final configuration of 338 
target objects, with a further 49 fibres assigned to guide 
stars and sky. 

The observations were scheduled on 2001 July 20, us- 
ing the 270R and 316R gratings. A total integration time 
of 1.5 hours was obtained in three 30- minute integrations 
to allow cosmic ray removal via a sigma-clipping algo- 
rithm. This gave a wavelength coverage of approximately 
3600-8000 A, varying slightly between the two spectro- 
graphs due to the slightly different dispersions of the grat- 
ings and from fibre to fibre on each spectrograph. The 
instrumental resolution was 9.9 A FWHM. 

4. Data Reduction 

Data were reduced using the 2df dr pipeline software pack- 
age, which was developed for the 2dF Galaxy redshift sur- 



vey (|Collessl[l99l IColless et al.ll200l[) . A full description 



of the semi-automated extraction and reduction procedure 
is given in these references. 

The software performs two independent redshift esti- 
mations depending on the emission and absorption proper- 
ties of the spectra. For the absorption spec tra we used the 
standard technique of cross-correlation ( "Tonrv fc Da^ 
1979i) with a set of eight absorption templates, five galax- 
ies and three stars. Th e galaxies used and th eir morpholog- 
ical types, as listed in IColless et al. 1 (l200lh . are: (1) NGC 
3379 (E), (2) NGC 4889 (cD), (3) NGC 5248 (Sbc), (4) 
NGC 2276 (Sc) and (5) NGC 4485 (Sm/Im). The stars 
used and their spectral classes are (6) HD 116608 (AlV), 
(7) HD 23524 (KOV) and (8) BD05°1668 (M5V). 

Spectra are prepared for cross-correlation via a six- 
step process that includes: continuum-subtraction; clip- 
ping strong emission lines; rebinning to a logarithmic 
scale; apodization; Fourier transformation; and applica- 
tion of an exponential fil ter. Furt h er de tails of each pro- 
cess can be found in .CoUess et al.l (|200l[) . 



The highest peak in the cross-correlation function is 
then fitted with a quadratic equation to obtain its position 
and height. The ratio of the peak height to the noise level 
in the cross-correlation function is subsequently computed 
and used to determine the quality of the redshift estimate, 
via the assignment of a quality fiag, Qa, to each spectrum. 
The value of Qa ranges from 1 to 4, with 4 and 3 meaning 
a reliable fit, 2 a probable fit and 1 an unreliable fit. If no 
fit is obtained a quality factor of is assigned. When Qa 
is 3 or 4, there is an additional requirement that at least 
four or six of the eight templates, respectively, produce the 
same redshift to within 600 km s~^. In addition, the code 
also outputs which of the eight templates produced the 
best correlation with each spectrum; this gives information 
on the object's morphological type. 



For emission spectra the code performs a Gaussian pro- 
file fit to each line after first subtracting the continuum. 
Any peaks in the resultant spectrum which are higher than 
3.3 times the RMS noise are marked as candidate emission 
lines. The three strongest lines in each spectrum are then 
tested in pairs for line separations consistent with a known 
redshift from a set of common emission lines ([Oil], H/3, 
[O III] , Ha or [N II] ) . If a match is found then all other 
emission lines that match that redshift to within 600 km 
are found, and the mean redshift is adopted as the 
emission redshift. If no two emission lines are found to 
match a particular redshift for a pair of common lines, (or 
there are fewer than three lines in total) then up to two 
single-line redshifts are retained for comparison with the 
absorption derived redshift. A quality flag for this redshift 
estimate, Qo, is assigned based on the number of features 
fitted. This number can take the values 0,1,2,4 with be- 
ing no lines, 1 for one line, 2 for two lines and 4 for three 
or more. 

The code then compares the redshift obtained via each 
method and selects the one with the higher quality flag. In 
the case where Qa = Qe the absorption value is selected. 
The best quality flag , Qb is then given as the greater of Qa 
or Qc, with two exceptions: (1) if the difference in redshifts 
derived via the two methods is less than 600 kms~^ then 
Qb is the greatest of Qa, Qo or 3 and (2) if Qa and Qo are 
both less than 2 and the difference in redshifts derived via 
the two methods is greater than 600 kms^^ then Qb is set 
to 1. 

After the automatic quality flags are assigned the spec- 
tra are then inspected manually and interactively re-fltted 
if necessary. Also at this time the observer may assign their 
own quality flag Q (independent of the automatic quality 
flag Qb) on a scale of 1 to 5, with 5 corresponding to 
spectra of the highest quality. 

5. Results 

Of the 338 spectra obtained, 136 turned out to be mis- 
classified stars. This was later tracked down to a problem 
with the object discriminator software for bright objects 
in the pre-release SuperCOSMOS data. For target selec- 
tion we set the lower magnitude limit to ttib — 15.5, at 
the bright end of the magnitude distribution of cluster 
members. Unfortunately, this was too bright for the ob- 
ject classification software to distinguish reliably between 
stars and galaxies. 

Figure [T] shows the distribution of redshifts obtained 
from all object fibres, with bin intervals of Az = 0.002. 
This figure has two main features; the misclassified stars 
near zero redshift and A3667 members at around z = 
0.056. The inset A3667 histogram shows remarkably lit- 
tle velocity structure either surrounding the main clus- 
ter group or within it. There is no evidence of distinct 
in-falling groups or other velocity features around the 
main cluster as found in other merging sy stems like the 



A3125/A3128 complex (jRose et al.l 120021 ). Then again. 



perhaps this is not surprising as clusters which harbour 
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A3667 Redshift Distribution Plot 



0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 0.22 

Redshift 

Fig. 1. 2dF redshift histogram for the A3667 field; the 
inset plot shows an enlarged version of the cluster redshift 
distribution. 



diffuse radio emission have been shown to be statistically 
more isolated than clusters of similar X-ray luminosity 
(jSchuecker fc Bohringer 1999). 

After the stars had been removed, the remaining 202 
spectra were examined individually for quality of fit. 
Spectra for which the quality flag Qb was less than 2 were 
treated with caution and 18 spectra with low signal-to- 
noise were removed from the sample. Of the remaining 
184 spectra, 143 fell within the range z = 0.044 - 0.068 
and are considered to be members of A3667. This redshift- 
restricted sample of 143 galaxies is published online at 
CDS in the format of Table [2l a sky plot of this sample, 
together with the 154 values from the literature, is shown 
in Figure [21 The remaining 41 galaxy spectra not related 
to A3667 but obtained during these observations are also 
published online at CDS in a format similar to those in 
Tabled 

There are 67 galaxies in common between the 2dF 
and literature samples. As can be seen from Figure O 
there is generally good redshift agreement between the 
two datasets. Analysis of the differences in redshift gives 
a median difference of zero and a standard deviation of 
Az = 0.00036 (102 kms"!), well below the combined red- 
shift uncertainty 200 kms""'^). This suggests that it is 
acceptable to combine the two datasets. 

In merging the two datasets we took note of the 2dF 
quality parameters. The blunder rate for spectra measured 
in the 2dFGRS for which 2df dr has assigned Qb > 2, is 
only 1.6% and this depends strongly on ma gnitude, with 



almo st all blunders occurring for me > 18 ijColless et al. 



20011 ). Thus, for a comparatively bright sample such as 
this, one can have confidence in the reliability of the 2dF 
redshifts. When combining the two datasets for subse- 
quent analysis the 2dF redshifts were used in preference to 
the literature values in all cases where Qb was greater than 
two. Of the nine 2dF spectra with Qb < 2, seven agreed 



3 -56.8 

O 
lU 

Q -57 




2Df data 
o Literature Data 



305.5 305 304.5 304 303.5 303 302.5 302 301.5 301 

RA (J2000) 

Fig. 2. Distribution on the sky of galaxies with measured 
redshifts in the region of A3667. Filled circles are the 2dF 
data and open circles are the redshifts from the literature. 
Where the 2dF and literature datasets agree in both po- 
sition and redshift, the circle is shown filled and outlined. 
The arcs denote the outside edges of the radio relics. 



well with the previously published redshifts and were pre- 
ferred for the combined dataset. The remaining two spec- 
tra showed large discrepancies of Az — 0.0069 and 0.0027 
(corresponding to ~ 2000 and 800 kms^^ respectively). 
These points can be seen on Figure [3] as outliers num- 
bered 105 and 78 (which correspond to their identifica- 
tion numbers in Table [21). Both objects were quite bright 
with magnitudes of 16.90 and 17. 39, and quoted e rrors 
of 61 and 200 kms~^ respectively ( Sodre et al.lll992h : for 
consistency, the 2dF values were retained for analysing 
the redshift distribution. The remaining two labelled out- 
liers, 70 and 101, were assigned the highest quality flags 
by 2df dr and were thus assumed to be correct. This gave 
a combined dataset of 231 redshifts in the range 0.044 to 
0.068, which were then used in the subsequent analysis. 

It should be noted that these observations combined 
with those previously published mean that A3667 has been 
spectroscopically observed at a completeness level of 70% 
down to our limiting magnitude of hj =18.3. Figure [4l 
shows all galaxies in a field of one degree radius centred 
on the peak in X-ray emission down to hj =18.3 as gray 
points overlaid with all galaxies for which spectroscopic 
redshifts are currently available. 



5.1. Redshift Distribution 

The redshift distribution for the combined dataset is 
shown in Figure [5l overlaid with a Gaussian fit to the 
central peak at z — 0.0555. While the redshift distribu- 
tion is well modelled by a Gaussian, there is some evi- 
dence for a small subgroup of galaxies at a slightly higher 
redshift; however, the galaxies belonging to this possible 
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Table 2. Redshifts for the first 10 galaxies* in the sample of 143 measured with 2dF in the cluster A3667 and falling 
in the redshift range 0.044 <z< 0.068 



ID 


RA 


Dec 


rriB 


2abs 




2avc 




morphological 




(J2000) 


(J2000) 












type 


1 


20 05 37.49 


--56 41 06.1 


16.89 


0.0558 








2 


2 


20 06 26.98 


-56 19 50.9 


18.16 


0.0586 








3 


3 


20 07 22.03 


-56 39 44.0 


16.41 


0.0535 








2 


4 


20 07 29.09 


-57 28 46.2 


18.09 


0.0553 


0.0553 


0.0553 




5 


5 


20 07 30.08 


-57 22 58.6 


16.75 


0.0521 


0.0518 


0.0520 




3 


6 


20 07 32.50 


-56 34 45.9 


17.43 


0.0560 


0.0557 


0.0559 




4 


7 


20 07 56.06 


-56 16 01.3 


16.65 


0.0543 


0.0545 


0.0544 


0.0535 


5 


8 


20 08 01.56 


-56 06 07.6 


17.60 


0.0594 


0.0594 


0.0594 




3 


9 


20 08 09.17 


-56 30 21.7 


16.80 


0.0523 






0.0524 


2 


10 


20 08 27.61 


-56 20 42.5 


17.79 


0.0572 


0.0569 


0.05705 
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The full data table contains 143 objects and is available in electronic form at CDS via anonymous ftp to cdsarc.u-strasbg.fr 
(130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A-|-A/(voI| /(page). Column 1 gives the source ID number; 
columns 2 and 3 are the Right Accession and Declination in J2000 coordinates; column 4 is the blue magnitude; columns 6 
& 7 are the redshifts as determined by absorption and emission features respectively and column 8 is the average redshift for 
cases where a redshift could be determined by both emission and absoprtion features; column 9 is the previously published 
redshift and column 10 is the morphol ogical type determin ed by the best correlation to the spectral templates used by 2dfdr, 
numbers correspond to those listed in lCoUess et al.l l|200ll ). 




0.043 0.048 0.053 0.058 0.063 

2dF Redshift 

Fig. 3. Redshift comparison for A3667. The literature red- 
shift is plotted against the 2dF redshift for the 67 sources 
in common. Clustering about the slope of unity suggests 
there is no overall systematic difference between the two 
datasets. The four labelled outliers are discussed in the 
text; the number corresponds to its identification number 
in the online version of Table [21 



subgroup are not spatially localised but are spread across 
the cluster. To search for possible subgroups a pseudo- 
3D plot of RA, Dec and redshift was constructed using 



-55.5 




-58 \ , , , , , , , , , 

305.5 305 304.5 304 303.5 303 302.5 302 301.5 301 

RA (J2000) (deg) 

Fig. 4. Coverage of galaxies for which spectroscopic red- 
shifts have been determined. Filled gray circles are all 
the galaxies down to hj < 18.3 in a one degree radius 
about A3667 as determined from SuperCOSMOS plate 
scans while open circles represent those galaxies for which 
redshifts have been obtained. The average completeness 
over the field is 70%. 

Matlab to explore the spatial distribution of all 231 galax- 
ies. While such a plot does not give a physical 3D model 
of the cluster, it is a useful tool in searching for associa- 
tions. Standard redshift versus ARA and ADec plots were 
also made (see Figure |6]) and tested for associations, on 
the premise that if subgroups are present there will be 
tight associations between the member galaxies in both 
plots. No likely subgroups were found via visual inspec- 
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0.052 0.056 

Redshift 



Fig. 5. Redshift distribution for A3667 with a Gaussian 
fitted to the central peak. The Gaussian is centred at the 
mean redshift of z = 0.0555 and there is a weak excess at 
around 0.0625. The velocity dispersion is 1102 km s~^. 



tion and the cluster appears extremely well mixed. This 
lack of velocity structure is consistent with previous re- 
sults (|Girardi et al.lll996l ). 

Statistical analysis of the total dataset gives a mean 
redshift of 0.0555, and a velocity dispersion, CTv, of 1102 
±46 kms~^, which is slightly higher than the previously 
published values. 



5.2. Quantitative Tests for Substructure 

In order to produce a more quantitative analysis 
of the likely substructure present w e performed 2 
and 3-dimensional Lee-Fitchett tests (|Fitchettl Il988t 
Fitchett fc: Websteii 119871 ) in addition to calculating the 
A statistic for the cluster ( Dressier fc ShectmanI 1988 ). 



5.2.1. Lee-Fitchett Test 



We performed a Lee-Fitc hett (LF) test (iFitchettI [l985 
Fitchett fc Websteiill987l ) on the 231 galaxies known to 
be members of A3667. The ratio obtained was 5.15 which 
is highly significant. However, as the LF test is simply 
a measure of spatial variations in variance, this does not 
imply subclustering in the case of A3667 which is highly 
aspherical. Although the Lee-Fitchett test has been shown 
to be a good indicator of substructure it does have limi- 
tations in that it will breakdown in the case of two widely 
spatially separated Gaussian velocity distributions. In the 
case of A3667 the LF test returns a result for significant 
substructure only in the plane of the sky along the axis 
of elongation of the galaxy positions. This is to be ex- 
pected given that the premise of the test is that clusters 
are roughly spherical. In order to assess the veracity of 
the test we simulated several clusters of the same size as 
A3667 with a high degree of aspherity (elliptically dis- 
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Fig. 6. Redshift distribution for the combined dataset of 
231 galaxies as a function of RA and Dec, measured in 
degrees from the cluster centre. 



tributed in RA-DEC) but which have a Gaussian velocity 
distribution. In all cases the simulated clusters also re- 
turned a highly significant substructure value for the LF 
test along the major axis of the ellipse defining the galaxy 
positions in RA & DEC. Neither in the case of A3667 
or the simulated clusters was there any significant axis of 
substructure associated with the velocity data. Thus, all 
we can say from this test is that distribution of galaxies in 
A3667 is elongated in the plane of the sky, as seen readily 
from the RA-DEC plots. 



5.2.2. A Statistic 

We undertook the calculation 



of the A Statistic 

( Dressier fc Shectman 19881) . which tests whether the dis- 
tribution of radial velocities {v, a) of a local {N = 15) 
population of galaxies changes with position across the 
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Table 3. KMM bimodal fitting results for the A3667 red- 
shift sample. 





Input Parameters 


Output Parameters 


Group 1 


Group2 


Group 1 


Group2 


Mean Group z 


0.0555 


0.0625 


0.0555 




Proportion 


0.935 


0.065 


1.0 





No. of Galaxies 


216 


15 


231 






cluster. The A value returned for A3667 was 1.26. For 
comparison 100,000 Monte-Carlo simulations, in which the 
positions of the galaxies were fixed in RA and Dec but ran- 
domly shuffled in radial velocity, returned a distribution 
of A values ranging from 0.76 to 1.49 with a mean of 1.09 
and standard deviation of 0.10. The lack of a clear dif- 
ference from the A value determined for A3667 suggests 
that there is no overall correlation between position and 
redshift in these data and, therefore, no significant sub- 
structure. 



5.3. Bimodality Testing 

A series of tests were performed to pursue the issue of 
bimodality in the cluster raised by P88 and S92. 



5.3.1. Velocity Gradient 

If the bimodal distribution reported by P88 is real and the 
two subgroups are at all inclined from the plane of the sky, 
we would expect to observe a gradient in radial velocity 
across the cluster. The Matlab cube referred to in Section 
15.11 was used to search for a velocity gradient across the 
cluster; none was found. 



5.3.3. Kurtosis and Skewness 

Tests for kurtosis and skewness were applied to the veloc- 
ity data. The kurtosis returned a value of 0.233 which is 
within 1 Standard Error of Kurtosis (SEK) for a Gaussian 
distribution with the same number of points (SEK = 
0.3223). The skewness of the velocity distribution was 
found to be 0.134 which again was within 1 Standard Error 
of Skewness (SES) for a Gaussian distribution with 231 
values (SES=0.16). For skewness to be considered signifi- 
cant we would expect values approaching two. Thus there 
appears to be no significant skewness or kurtosis for the 
velocity distribution as compared with a Gaussian distri- 
bution of the same number of data points. 

5.3.4. Isodensity Plots 

In order to assess the bimodality of the galaxy surface 
distribution, Gaussian smoothing was applied to the spa- 
tial distribution of the 231 galaxies with redshifts be- 
tween 0.044 and 0.068. Figure [7] shows the galaxy loca- 
tion data smoothed with Gaussian filters ranging in size 
from cr = 0.20° to 0.35°. The 0.20° filter image shows an 
excellent overall match to the published isodensity plot of 
P88, except that the rcdshift-restricted data do not exhibit 
the significant subgroup structure seen by P88. However, 
there is a definite elongation of the cluster and some ev- 
idence of a slight enhancement in galaxy density in the 
NW region where P88 find a significant subgroup. There 
is also a feature of equal significance on the SE side of 
A3667. This second feature is in the approximate loca- 
tion that Ivikhlinin fc MarkevitchI (j2002h predict a small 
subgroup based on their Chandra analysis of the X-ray 
features. However, it is likely that both these features re- 
sult from the limited sample size rather than comprising 
distinct subgroups. Neither feature persists if the filter size 
is increased even marginally. 



5.3.2. Kaye's Mixture Model 

To test quantitatively whether a single Gaussian is an ad- 
equate fit to the redshift distribution shown in Figure 
El we used Kaye's Mixture Model algorithm, KMM 



([Ashman et al.lll994l ). which allows a user-specified num- 
ber of Gaussian profiles to be fitted to the unbinned, raw 
data and provides a statistical comparison with a single 
Gaussian fit. The user provides an estimate of the mean 
value of each peak and the fraction of the members in each 
peak. The software, provided by Keith Ashman, then as- 
signs membership of each redshift to one of the initially 
supplied Gaussian profiles, and returns posterior proba- 
bilities of group membership for each profile. 

The algorithm was applied for the two-profile case and 
Table [3] summarises the input and output parameters. For 
the 231 redshifts measured for A3667, the code assigned 
all values to be members of the main group centred about 
z = 0.0555, with an overall confidence level of 96.7%. 



5.4. Galaxy Morphologies 

As mentioned above, one of the products of the 2d.fd.r 
cross-correlation redshift determination is the template 
which gives the strongest correlation with each measured 
spectrum. This provides information on each galaxy's 
morphological type. Examination of the spatial distribu- 
tion of various galaxies within the cluster may provide 
valuable information on the cluster's dynamical history. 

Figure [8] shows the distribution of galaxy types in 
A3667 as defined by the 2df dr template. The morpholog- 
ical types are elliptical (E), central dominant (cD), spiral 
(Sc, Sbc) and irregular/small (IM/SM). Not surprisingly 
the predominant galaxy type is elliptical. The spiral popu- 
lation occurs mainly in the north-western part of the clus- 
ter and at first glance this appears to be significan t. We ap- 
plied a bi-dimensional Kolmogorov-Smirnoff test ([Peacock 
19831 ) to test if there was a statistically significant dif- 
ference between the distribution of spiral and elliptical 
galaxies in the cluster. Using positions of the 85 elliptical 
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Fig. 7. Galaxy isodensity plots for A3667, generated by 
application of Gaussian filters of various sizes (from a = 
0.2 to 0.35°) to the 231 known cluster members. Contours 
are shown at V2 intervals. 

galaxies and 19 spirals with morphological types identi- 
fied in the 2dF observations we computed a P(> Zoo) of 
0.999 which supports the null hypothesis that the distri- 
butions are drawn from the same population. In addition 
we calculated the mean redshift and velocity dispersion 
for the separate populations finding them to be almost 
identical with Zdiipticai = 0.05570 and ^spiral = 0.05572 
while (7cnipticai=1007 ± 109 kms~^ and aspirai=1133 ± 
260kms-i. 
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Fig. 8. Galaxy type spatial distribution for the 143 
ies in the 2dF dataset. 
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Observations of starburst galaxies have shown a strong 
correlation between star formation and tidal disruption. 



Such disruption may be either the result of galaxy merg- 
ers, such as the famous case of the Antennae galaxies 
(NGC4038 & NGC4039) (jWhitmore et al.lll999h . or from 
accretio n or merger shoc ks in the intra-cluster medium 
(ICM) ( Burns et al.ll200^ . An excess of starburst galax- 
ies would seem then to indicate a turbulent state, either 
through more in dividual galaxy coUisfons a s in the case 



of CL 1358-h62 ([Bartholomew et al.l l200lh . or through 



Shockwave disruption throughout the ICM, possibly re- 
sulting in ram pressure stripping. Thus, one might ex- 
pect that an overdensity of starburst galaxies would be 
evidence for a recent cluster merger. However, during 
the late stages of cluster evolution we know from the 
Butcher-Oemler effect that cluster constituent galaxies are 
more likely to be gas-stripped red ellipticals than gas- 
rich blue spirals. Without an abundance of gas it is un- 
likely then that starbursts per se would be observed. A 
more likely indicator of recent merging might be an ex- 
cess of the so-called 'E-l- A' galaxies, which repres ent a 
post-starburst population ( Couch fc Sharpie^ 1987 ) gen- 
erated either throug h galaxy-galax y mergers or ram pres- 
sure stripping (jBelloni et al. I [1991). Several studies have 
addressed the prevale nce of post-starburst galaxies iii 
cluster enyironin ents ( Hogg et al. 200d_ Pogg ianti 20031: 
Tran et al.l l2007f ) as weh as the field (|Blake et al. 200l 



finding that the fraction of post-starburst galaxies is dif- 
ferent to the field at low redshifts (z< 0.1). As an excess 
of post-starburst galaxies may be an indicator of quench- 
ing of star for mation in galax ies as a result of the cluster 
environment ( Poggiantil 20031 ) they can be considered an 

important indicator or merger activity. 

No E+A galaxies (|Couch fc Sharplesll 19871 ) were found 
among the spectra. Perhaps this is not surprising as only 
37 examples with strong H/3, H7 and H(5 lines were found 
in a subset of over 15,000 low- to mid-redshift galaxies 
observed in the 2dFGRS. This is an incidence rate of 
only 0.25 percent (Deeley 2000, private communication). 
While at moderate to high redshifts E-t-A galaxies would 
be more likely to be found in regions where galaxy-galaxy 
i nteractions are p rominent, i.e., the central cluster region 
( Tran et al.ll2003l ) at redshifts less tha n one they tend to be 



foimd in the field ( Blake et al. 2004[ l an d only a few have 
to date been seen in low redshift clusters (jCaldwell fc Rose 
1997t ICaladbOOOl ) like A3667. Due to the spacing hmita- 
tions imposed by the physical size of the 2dF fibre buttons 
the central region of A3667 is undersampled by 2dF, as 
can be seen in Figure [2] meaning that if there were E-l-A 
galaxies in the cluster we would potentially miss them. 

6. Discussion & Conclusions 

Cluster mergers are rarely established based solely on ev- 
idence from one wavelength regime, generally requiring a 
multiwavelength approach to reveal the true dynamical 
nature of the interactions. 

The presence of diffuse radio emission in clusters 
of galaxies points to the existence of cluster-wide mag- 
netic fields of the order of 0.1 -1 /xG and a population 
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of relativistic electrons with Lorcntz factor, 7 ^ 1000 
and e nergy densities of IQ-^^ - 10"" erg cm ([Ferettil 
2004 ^ . Diffuse radio emission is observed most commonly 
in clusters which are believed to have undergone a re- 
cent merger event and this leads naturally to the con- 
clusion that the phenomenon is due to turbulence and 
shocks in th e intra-cluster medium energizing the electron 



Several authors h ave claimed that A3667 is observed in 
a pos t- merger state ( Knopp et alj 19961: Markevitch et al 



snocks m th e mtra-ciuster medium energizing tne electron 
population (^Giovannini fc Feretti|[2003iHarris et al.lll980 : 



Tribblelll99.C Simulations of the turbulence and shocks in 
the cluster environment are able to pre dict the characteris- 
tics of observed diffuse radio emission ( [ Miniati et al]|2001 



Pfrommer et al.l2007 ; Roettiger et al.lll999l ) and studies of 
the number counts suggest that the most likely lifetimes of 



diffuse radio emission is of the order of 1 Gyr (|Kuo et al 



20041) . Additionally, theoretical models of merging sys- 
tems suggest that X-ray substructure will be erased af- 
ter the sound crossing time and thus, clusters with dis- 
turbed X-ray isophotes are likely to be dynamically young 
and in on-going merger events. Finally, the strongest ev- 
idence for clusters being in recent dynamical interactions 
is the detection of temperature chang es across the source 
(|Govoni et al.ll2004l : IPratt et al.ll2007l i. 



While the number of clusters containing d iffuse emis- 
sion is still comparatively low ( Ferettil 2004), there are 
now a number of well studied systems with mul tiwave- 
length data inclu ding A2256 (IC arke fc Ensslinl [ioOQl . 
A754 ()Henrv et al.ll2004l) . IC1262 (iHudson et al.ll2o61as 
well a. s the numerous examples presented by Govoni et al.l 
(j2004r ). To date most analysis tends to concentrate on the 
radio and X-ray properties of the clusters. However, often 
this is insufficient to determi ne the merger geometry. For 
example, the case of A3921 ( Ferrari et ahllioOQ) required 
data from optical, radio and X-ray analysis to determine 
the dynamical state of the cluster. In this instance optical 
observations alone were inconclusive but taken together 
with the other wavelength data the merger axis could be 
established. Similarly, high resolution X-ray from XMM- 
Newton analyzed in conjunction with existing multiwave- 
length data has proved vital to unrav eling the cluster dy- 



namics in s everal cases such a s A85 (jPurret et al 
and A3921 (|Belsole et al.ll2005[ ). 



2005h 



A3667 has a mean redshift of 0.0555, a high velocity 
dispersion (1102 ± 46 kms^^) and an elongated optical 
axis. It is typical of rich clusters in that it is dominated by 
elliptical galaxies. There is marginal evidence for a group- 
ing of spirals in the north-western part of the cluster. The 
velocity distribution is well modelled by a single Gaussian 
peaking at z = 0.0555, with a weak excess around a red- 
shift of 0.0625. There seems to be little or no support for 
distinct in-falling groups either from the velocity informa- 
tion or the spatial distribution of galaxies with measured 
redshifts. Moreover, there is no evidence for a velocity gra- 
dient across the cluster suggested by the previously re- 
ported bimodal optical distribution. However, A3667 is a 
rich cluster and to date only 231 redshift have been ob- 
tained for cluster members. A larger sample is needed to 
categorically rule out the presence of substructure. 



ig^ 

iVikhHnin et al.ll200ll ). This may be consistent with the 
high velocity dispersion and mixing in the cluster. If 
the observed diffuse radio emission was generated via a 
cluster-cluster merger it is expected that the collisional 
velocities would be a few thousand kms^"'^. If such a 
merger occurred, the current data suggest that the merger 
axis must be very close to the plane of the sky. This 
would partly account for the lack of observable subgroup 
peaks in the velocity data. 

Full understanding of the the dynamical situation in 
A3667 will require multi- wavelength analysis. We have ob- 
tained detailed new radio data for the entire cluster at 1.4, 
& 2.4 GHz, we will use this in conjunction with previously 
published radio (843 MHz) and X-ray data to further in- 
vestigate this cluster. Findings will be published in a fu- 
ture paper. 



Il99 9t lRoettig;er et all Il999t iRottgering et all Il997 
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